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Abstract 

The problems encountered in trying to quantize the various cosmological mod- 
els, are brought forward by means of a concrete example. The Automorphism 
groups are revealed as the key element through which G.C.T.'s can be used for 
a general treatment of these problems. At the classical level, the time dependent 
automorphisms lead to significant simplifications of the line element for the generic 
spatially homogeneous geometry, without loss of generality. At the quantum level, 
the "frozen" automorphisms entail an important reduction of the configuration 
space -spanned by the 6 components of the scale factor matrix- on which the 
Wheeler-DeWitt equation, is to be based. 

In this spirit the canonical quantization of the most general minisuperspace ac- 
tions -i.e. with all six scale factor as well as the lapse function and the shift vector 
present- describing the vacuum type II, I geometries, is considered. The reduction 
to the corresponding physical degrees of freedom is achieved through the usage of 
the linear constraints as well as the quantum version of the entire set of all classical 
integrals of motion. 
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1 Introduction 



Since the conception by Einstein of General Relativity Theory, a great many efforts have 
been devoted by many scientists to the construction of a consistent quantum theory of 
gravity. These efforts can de divided into two main approaches: 

(a) perturbative, in which one splits the metric into a background (kinematical) part 
and a dynamical one: g^ u = rj^+h^ and tries to quantize h^. The only conclusive 
results existing, are that the theory thus obtained is highly nonrenormalizable |lj]. 

(b) non perturbative, in which one tries to keep the twofold role of the metric (kinemat- 
ical and dynamical) intact. A hallmark in this direction is canonical quantization. 

In trying to implement this scheme for gravity, one faces the problem of quantizing a 
constrained system. The main steps one has to follow are: 

(i) define the basic operators "g^ and 7r M1/ and the canonical commutation relation 
they satisfy. 

(ii) define quantum operators whose classical counterparts are the constraint func- 
tions H^. 

(iii) define the quantum states ^[g] as the common null eigenvector of H^, i.e. these 
satisfying H^[g] = 0. (As a consequence, one has to check that H^, form a closed 
algebra under the basic CCR.) 

(iv) find the states and define the inner product in the space of these states. 

It is fair to say that the full program has not yet been carried out, although partial steps 
have been made @. 

In the absence of a full solution to the problem, people have turned to what is 
generally known as quantum cosmology. This is an approximation to quantum gravity 
in which one freezes out all but a finite number of degrees of freedom, and quantizes the 
rest. In this way one is left with a much more manageable problem that is essentially 
quantum mechanics with constraints. Over the years, many models have appeared in the 
literature ||. In most of them, the minisuperspace is flat and the gravitational field is 
represented by no more than three degrees of freedom (generically the three scale factors 
of some anisotropic Bianchi Type model Q). 

In order for the article to be as self consistent as possible, we include in section 2, a 
short introduction to the theory of constrained systems and in section 3, the Kantowski- 
Sachs model is treated both at the classical and the quantum level -as an interdisciplinary 
example. In section 4, the importance of the Automorphism group is uncovered and the 
quantization of the most general Type II, I Vacuum Bianchi Cosmologies, is exhibited. 
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2 Elements of Constrained Dynamics 



2.1 Introduction 

In these short notes, we present the elements of the general methods and some techniques 
of the Constrained Dynamics. It is about a powerful mathematical theory (a method, 
more or less) -primarily developed by P. A. M. Dirac. The scope of it, is to describe 
singular (the definition is to be presented at the next section) physical systems, using a 
generalization of the Hamiltonian or the Lagrangian formalism. This theory, is applicable 
both for discrete (i.e. finite degrees of freedom) and continua (i.e. infinite degrees of 
freedom) systems. 

For the sake of simplicity, the Hamiltonian point view of a physical system is adopted, 
and the discussion will be restricted on discrete systems. A basic bibliography, at which 
the interested reader is strongly suggested to consult, is quoted at the end of these notes. 
Also, the treatment follows reference ||. 

2.2 The Hamiltonian Approach 

Suppose a discrete physical system, whose action integral is: 



The dynamical coordinates, are denoted by q l , with i £ [1, . . . ,N] The Lagrangian is a 
function of the coordinates and the velocities, i.e. L = L(q\q l ). 

A note is pertinent at this point. If one demands the action integral Q2.2.1|) , to be 
scalar under general coordinate transformations (G.C.T.), then he can be sure that 
the content of the theory to be deduced, will be relativistically covariant even though 
the form of the deduced equations will not be manifestly covariant, on account of the 
appearance of one particular time in a dominant place in the theory (i.e. the time 
variable t occurring already, as soon as one introduces the generalized velocities, in 
consequently the Lagrangian, and finally the Lagrange transformation, in order to pass 
from the Lagrangian, to the Hamiltonian). 

Variation of the action integral, gives the Euler-Lagrange equations of motion: 



In order to go over to the Hamiltonian formalism, the momentum variables pi, are 
introduced through: 



In the usual dynamical theories, a very restricting assumption is made; that all momenta 
are independent functions of the velocities, or -in view of the inverse map theorem for a 






(2.2.2) 




(2.2.3) 
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function of many variables- that the following (Hessian) determinant: 

f) 2 T 

\H U \ = I -^-A I (2.2.4) 

is not zero in the whole domain of its definition. If this is the case, then the theorem 
guarantees the validity of the assumption, permits to use the Legendre transformation, 
and the corresponding physical system is called Regular. If this is not the case (i.e. some 
momenta, are not independent functions od the velocities), then there must exist some 
(say M) independent relations of the type: 

<pm(q,p) = 0, me[l,...,M) (2.2.5) 

which are called Primary Constraints. The corresponding physical systems, are charac- 
terized as Singular. 

Variation of the quantity piq 1 — L (the Einstein's summation convention is in use), 
results in: 

5 { Pl q l — L) — ... — (5 Pi ) q 1 - <V (2-2.6) 

by virtue of ( ^2~2|) . One can see that this variation, involves variations of the g's and 
the p's. So, the quantity under discussion does not involve variation of the velocities and 
thus can be expressed in terms of the g's and the p's, only. This is the Hamiltonian. It 
must be laid stress on the fact that the variations, must respect the restrictions ( |2.2.5| ), 
i.e. to preserve them -if they are considered as conditions (see, e.g. C. Caratheodory, 
"Calculus of Variations and Partial Differential Equations of the First Order", AMS 
Chelsea (1989)). 

Obviously, the Hamiltonian is not uniquely determined for, zero quantities can be 
added to it. This means that the following: 

H T = H + u m <p m (2.2.7) 

where w m 's are arbitrary coefficients in the phase space (including the time variable), is 
a valid Hamiltonian too. Variation of ( |2.2.7| ) results in: 



■i = an + md± + t th t van i s hes as (2.2.5) 

^ dpi dpi ' v i (2.2.8) 

Pi = ~ 1T~ ~ um ~d ~ term that vanishes as ( |2.2.5|) 



These are the Hamiltonian equations of motion for the system under consideration. This 
scheme, reflects the previous observation about variations under which, conditions must 
be preserved. 

In order to proceed, a generalization of the Poisson Brackets must be introduced. 
This is done as follows: 

Let f,g,hbe quantities on a space, endowed with a linear map { , } such that: 

{f, 9} + {9, /} = Antisymmetry 

{/ + 9,h} = {f, h} + {g, h} Linearity , , 

{f 9l h} = f{g, h} + {/, h}g Product Law 1 ' 

{/, {9, h}} + {9, {h, /}} + {h, {/, g}} = Jacobi Identity 
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If the space is the phase space, then these Generalized Poisson Brackets, reduce to the 
usual ones: 



{f,g} 



df dg 
dq 1 dpi 



dg df 
dq 1 dpi 



(2.2.10) 



otherwise are subject to the previous algebra -only. 

For a dynamical variable -say g, one can find -with the usage of: 

■ _ dg { dg . 

dq 1 dpi 1 

and of Q2.2.8| ), as well as the generalized Poisson Bracket Algebra (2.2.9): 

g w {g, H T } 



(2.2.11) 



(2.2.12) 



The symbol ~ is the Weak Equality symbol and stands for the following rule (deduced 
from a thorough analysis of the previous procedure): 

A constraint, must not be used before all the Generalized Poisson Brackets, are calculated 
formally (i.e. only with the usage of the algebra (\2. 2. 1 ) and the usual definition ( \2.2.1(J[) 
-when the last is applicable). This rule, is encoded as: 



m (g,p)^O, m£[l,...,M] 



(2.2.13) 



In the previous procedure, the position of that rule, reflects the need to manipulate the 
w m 's, which may depend on t only -since they are unknown coefficients, the definition 
( p.2. 10| ) can not be used. 

If the dynamical variable g is any one of the constraints, then (|2.2.5|) declare the 
preservation of zero. Thus, consistency conditions, are deduced: 

{0 m ^#} + M m {</w,<M = O (2.2.14) 

There are three possibilities: 

GC\ Relations ( |2.2.14j ) lead to identities -maybe, with the help of ( |2.2.5|) . 

CC<i Relations ( |2.2.14j ) lead to equations independent of the w's. These must also be 
regarded as constraints. They are called Secondary, but must be treated on the 
same footing as the primary ones. 

CC3 Relations ( |2.2.14| ) impose conditions on the u's. 

The above procedure must be applied to all secondary constraints. Again, the possible 
cases will be the previous three. The new constraints which may turn up are called 
secondary too. The procedure is applied for once more and so on. At the end, one 
will have a number of constraints (primary plus secondary) -say J- and a number of 
conditions on the it's. A detailed analysis of the set of these conditions, shows that: 



u m = U m (q,p)+V a (t)V a m (q,p) 



(2.2.15) 
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where V^ m (g,p) are the a (in number) independent solutions of the homogeneous systems: 

K m (g,p){0m»,<M = o 

The functions V a {t) are related to the gauge freedom of the physical system. 

Some terminology is needed at this point. 
A dynamical variable R, is said to be First Class, if it has zero Poisson Bracket, with all 
the constraints: 

{R,0 n } = O, ne[l,...,J] (2.2.16) 

where J is the total number of constraints -i.e. primary plus all the secondary ones. 
It is sufficient for these conditions, to hold weakly -since, by definition, the 0's are the 
only independent quantities that vanish weakly. Otherwise, the variable R, is said to be 
Second Class. If R is First Class, then the quantity {R, <fr n } is strongly equal to some 
linear combination of the </>'s. The following relative theorem (with a trivial proof) holds: 
"The Poisson Bracket of two First Class quantities, is also First Class". 

Using the result ( |2.2.15| ) the Hamiltonian ( |2.2.7| ) , which is called Total Hamiltonian, 
is written: 

H T = H + U m <j) m + V a V™$ m = H' + V^a (2.2.17) 

with obvious associations. It can be proved that H' and (f) a , are first class quantities. 
With this splitting and the relation ( [2.2.12|) for a dynamical variable g, it can be deduced 
that: 

The First Class Primary Constraints 4> a , are the generating functions (i.e. the quantities 
{g,(p a }) of infinitesimal Contact Transformations; i.e. of transformations which lead to 
changes in the q 's and the p 's which do not affect the physical state of the system. 

Successive application of two contact transformations generated by two given First 
Class Primary Constraints and taking into account the order, leads -for the sake of 
consistency- to a new generating function: {g, {<f) a , 4> a '}} • Thus one can see that First 
Class Secondary Constraints, which may turn up from {0 a ,0 a /}, can also serve as gen- 
erating functions of infinitesimal Contact Transformations. Possibly, another way to 
produce First Class Secondary Constraints, is the First Class quantity {H', <p a }. Since 
no one has found an example of a First Class Secondary Constraint, which affects the 
physical state when used as generating function, the conclusion is that all First Class 
quantities, are generating functions of infinitesimal Contact Transformations. Thus, the 
total Hamiltonian should be replaced by the Extended Hamiltonian He, defined as: 

H E = H T + U a "<p a „ (2.2.18) 

where the a "'s are those First Class Secondary Constraints, which are not already 
included in H?. Finally, the equation of motion for a dynamical variable g ( ^.2.12| ) is 
altered: 

9~{9:H E } (2.2.19) 
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2.3 Quantization of Constrained Systems 
2.3.1 No Second Class Constraints are Present 

The quantization of a classical physical system, whose Lagrangian, gives first class con- 
straints only, is made in three steps: 

51 The dynamical coordinates g's and momenta p's, are turned into Hermitian Oper- 
ators q 's and p 's, satisfying the basic commutative algebra: = idj. 

5 2 A kind of a Schrodinger equation, is set up. 

53 Any dynamical function, become Hermitian Operator -provided that the ordering 
problem is somehow solved. 

Obviously, the constraints -being functions on the phase space- are subject to the S3 
rule. Dirac, proposed that when the constraints are turned into operators, they must 
annihilate the wave function \1/: 

Ji* = 0, V % (2.2.20) 

Successive application of two such given conditions and taking into account the order, 
for sake of consistency, results in: 

= (2.2.21) 

In order for operational conditions (|2.2.21|) not to give new ones on one demands: 

[Mj] =C*.,J fc (2.2.22) 

If it is possible for such an algebra to be deduced, then no new operational conditions on 
\1/ are found and the system is consistent. If this is not the case, the new conditions must 
be taken into account and along with the initial ones, must give closed algebra, otherwise 
the procedure must be continued until a closed algebra is found. The discussion does not 
end here. Consistency between the operational conditions ( |2.2.20| ) and the Schrodinger 
equation, is pertinent as well. This lead to: 

[Ji,H]* = (2.2.23) 

and consistency know, reads: 

(2.2.24) 
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2.3.2 Second Class Constraints are Present 

Suppose we have a classical physical system, whose Lagrangian, gives second class con- 
straints. Any set of constraints, can be replaced by a corresponding set of independent 
linear combinations of them. It is thus, in principle, possible to make arrangement such 
that the final set of constraints, contains as much first class constraints as possible. Us- 
ing the remaining -say S in number- second class constraints, the following matrix is 
defined: 

*ij = \Xi,Xj], (hj) € [1,... ,S] (2.2.25) 

where , are the remaining (in classical form) second class constraints. A theorem can 
be proved : 

"The determinant of this matrix does not vanish, not even weakly". 

Since the determinant of A is non zero, there is the inverse of this matrix; say A -1 . 
Dirac, proposed a new kind of Poisson Bracket, the { , }d- 

s s 

= Xi>Ay > (2-2.26) 

i=l j=l 

These Brackets, are antisymmetric, linear in their arguments, obey the product law and 
the Jacobi identity. It holds that: 

{g,H E } D n{g,H E } (2.2.27) 

because terms like {x^He}, with H E being first class, vanish weakly. Thus: 

gn{g,H E } D (2.2.28) 

But: 

{^ Xs } D = ... = (2.2.29) 

if £ is any of the g's or the p's. Thus, at the classical level, one may put the second class 
constraints equal to zero, before calculating the new Poisson Brackets. That means that: 

Mi The equations X — may be considered as strong equations. 

M 2 One, must ignore the corresponding degrees of freedom and 

M 3 quantize the rest, according to the general rules, given in the previous section. 
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3 A Pedagogical Example: 

The Kantowski-Sachs Model 



The purpose of the present section is twofold: 

• to illustrate -at the classical level- an application of Dirac's method for constrained 
systems. 

• to present, in an easy manner, the problems rased by the quantization of such 
systems. 

The example chosen, is that of Kantowski-Sachs reduced Lagrangian -i.e. of a vacuum 
cosmological model; thus the interdisciplinary character of the section, emerges. 



3.1 The Classical Case 

Consider, the Kantowski-Sachs model (described in ||), characterized by the line ele- 
ment: 

ds 2 = -N 2 (t)dt 2 + a 2 (t)dr 2 + b 2 (t)d6 2 + b 2 (t)stn 2 (6)d(f) 2 (3.3.1) 
The corresponding Eintein's Field Equations, are: 



Go °-~ \~m) ~ 2 imr~\m ] (3 - 3 - 2) 

i(f) \ 2 ( a(t)b'(t) \ 2 _ a(t)*V(t)N'(t) 2 a(t) 2 b"(t) 



_ b(t)a'(t)b'(t) a'(t)b(t) 2 N'(t) b(t)b'(t)N'(t) b(t) 2 a"(t) b(t)b"(t) 
22 a(t)N(t) 2 N(t) 3 a(t) N(t) 3 a(t)N(t) 2 N(t) 2 { ' ' ] 



G 33 = sin{6) 2 G 22 (3.3.5) 

The first of these (Goo); is the quadratic constraint equation i.e. its time derivative 
vanishes -by virtue of the other two (Gn,G'22)- This is a "peculiarity" of Einstein's 
system, and reflects the time reparametrization invariance t — > t — f(t). Under such a 
transformation, a(t) and b(t) change as scalars (a(t) = a(t), ditto the b(t)) while N(t), 
changes as density (N(t)dt = N(t)dt), revealing its nature, as a Lagrange multiplier. 

It must be brought to the reader's notice that the above set of equations (G^), can 
be obtained from the following action principle: 

S = I Ldt = I | _»(«)*■(«) + jg^M'l j (3.3.6) 
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-where a(t), b(t) and N(t), are the three degrees of freedom (the g's)- which has the 
above mentioned reparametrization invariance. 
The momenta are: 

v ~— ~ mHt) (3 3 7a) 

Pa ~dd(t) ~ N(t) (3 ' 3 - 7a) 

dL a(t)b(t) + a(t)b(t) 

Pb = W) = W) — (3 - 3 - Tb) 

dL 

p N = —. = 3.3.7c 

v dN{t) y ' 

From the third of (|3.3.7a|) , one can see that there is one primary constraint: 

(3.3.8) 

The total Hamiltonian is: 

H T = H + u(t)p N (3.3.9) 

where: 

H = p a a(t) + p b b(t) — L = N(t)tt(t) (3.3.10) 

with: 

The consistency condition ( |2.2.14| ) applied to: 
A\ the constraint p^ ~ 0, gives one secondary constraint: 

X = {p N , H} = {p N , N(t)tt(t)} = -n(t) « (3.3.12) 
A straightforward calculation, results in: 

{ X ,Pn} = (3.3.13) 

A 2 the previously deduced secondary constraint x ~ 0, gives -by virtue of ( |3.3.11| ), 
( |3.3.12|) and (|3.3.13| )- no further constraints, since it is identically satisfied (CC± 



case) : 

{ X ,H} + u(t){x,p N } = (3.3.14) 
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The Poisson Bracket ( |3.3.13| ) also declares that both and Xi are fi rs t class quan- 
tities. 

Finally, the equations of motion are: 



a(t) « {a(t),H T } 


(3.3.15) 


Pa ~ {Pa, H T } 


(3.3.16) 


b\t)*{b{t),H T } 


(3.3.17) 


p b ta {p b , H T } 


(3.3.18) 


N(t)*{N(t),H T } 


(3.3.19) 


Pn ~ {pn, H t } 


(3.3.20) 


The first four equations constitute the usual set of the Euler-Lagrange equations for the 
a(t) and b(t), degrees of freedom. Equation (|3.3.19|), results in the gauge freedom related 



to N(t) since -according to this equation- N(t) = u(t), i.e. an arbitrary function of 
time, while equation ( |3.3.20|) is trivially satisfied, in view of (|3.3.12|) . 



Finally a remark concerning the existence of shift terms of the form Ni(x 3 ,i)dx l dt 
-x 3 stand for (r,6,<p)- where Ni(x,t) = N a (t)af(x): their existence entails constraint 
equations (Go*) -again preserved in time, by virtue of the (Gy) equations- which reflect 
the space reparametrization invariance x l — > x l {x\ t). Along with the existence of these 
shift terms, a change in the spatial part of the line element, is induced. 

3.2 The Quantum Case 

In trying to quantize the previously described constraint Hamiltonian system, various 
problems, arise []5], |5|. 

In the canonical approach [0] -and references therein-, the Schrodinger representa- 
tion, is most frequently adopted. Applied to our example, this entails the step: 

a — > a = a 
b^b = b 
N -> N = N 

Pa^Pa = ~i 

Pb^Pb = -i 



Pn ^Pn = -i^j 



(3.3.21) 

db 



When trying to implement Dirac's proposal (steps Si, 5*2 of the section 2.3.1) we 
came across the factor ordering problem (see e.g. T. Christodoulakis, J. Zanelli, Nuovo 
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Cimento B 93 (1986) 1). Its resolution is achieved via the recipe that the kinetic term 
must be realized as the conformal Laplacian. This is due to the covariance in the change 
N(t) = N (t) f (a(t) , b(t)) -with the understanding that / (a(t) , b(t)) , is identified to an 
arbitrary function of time. The conformal Laplacian must be based on the metric: 

„ _ / a/26 2 -1/26 \ 
9 ~ V -1/26 J 

because of the correspondence principle, since H = N^g^piPj + V), where pi = p a , 
p 2 = Pb and V = —a/2. In two dimensions, the conformal Laplacian, reduces to the 
typical one (see next section for details). Thus, following Dirac's quantization program, 
we deduce: 

H T V = (3.3.23) 



(3.3.22) 



or: 



p N V = (3.3.24) 

(Constraint) and: 

= (3.3.25) 
(Wheeler-DeWitt equation). Under the transformation: 



a — > u = 6 
6 — > v — a 3 b 

the Wheeler-DeWitt equation, assumes the form 

. d 2 m 



dudv 

and under a second transformation: 



the Wheeler-DeWitt equation, takes the form: 

d 2 m d 2 ^ 



OX 2 dY 2 



(3.3.26) 



* = (3.3.27) 



u(t) ^X = ^ 
v (t) ^ Y = ^ 



- * = (3.3.29) 



Now the previous equation can be solved via the method of separation of variables, e.g. 
^(X, Y) = A(X)B(Y); its general solution, consists of products of Exponentials and/or 
Trigonometric functions -depending on the sign and the value of the separation constant. 

Of course, in order to complete the program of quantization, we need to construct 
the Hilbert space, i.e. to select a measure. The problem is open, because there is an 
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infinitude of candidates. If one invokes some sort of " naturality" , one could adopt as 
a measure the square root of the determinant of the supermetric, i.e. 2b -in our case. 
This however, causes two unpleasant drawbacks: the first is that the wave function, is 
not square integrable, and the second is the violation of the conformal covariance. 

In the case where shift terms and more spatial metric cross terms, are present, one 
would like to know, what features of the above exhibition, are generic and thus, char- 
acterize the general situation. The answer is given through the consideration of the au- 
tomorphism group, which can be considered as the symmetries of the symmetry group 
of the 3-space. Their action entails considerable simplification, both at the classical and 
the quantum level. The spirit of these ideas, is exhibited in the next sections. 
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4 Automorphisms in Classical & Quantum Cosmology 



4.1 The Simplification of Einstein's Equations 

It has long been suspected and/or known, that automorphisms, ought to play an impor- 
tant role in a unified treatment of this problem. The first mention, goes back to the first 
of |J. More recently, Jantzen, -second of §J- has used Time Dependent Automorphism 
Matrices, as a convenient parametrization of a general positive definite 3x3 scale fac- 
tor matrix r y a pif)i in terms of a -desired- diagonal matrix. His approach, is based on 
the orthonormal frame bundle formalism, and the main conclusion is (third of M, pp. 
1138): "... the special automorphism matrix group SAut(G), is the symmetry group of 
the ordinary differential equations, satisfied by the metric matrix 7^, when no sources 



are present ..." Later on, Samuel and Ashtekar in ||10|| , have seen automorphisms, as 
a result of general coordinate transformations. Their spacetime point of view, has led 
them, to consider the -so called- " Homogeneity Preserving Diffeomorphisms" , and link 
them, to topological considerations. 



4.1.1 Time Dependent Automorphism 
Inducing Diffeomorphisms 

It is well known that the vacuum Einstein field equations can be derived from an action 
principle: 

A = ^- I V- (i) 9 {i) R d 4 x (4.4.1) 

l07T J 

(we use geometrized units i.e. G — c — 1) 

The standard canonical formalism |TT] makes use of the lapse and shift functions ap- 
pearing in the 4-metric: 

ds 2 = (N*Ni - N 2 )dt 2 + 2N i dx i dt + g ij dx i dx j (4.4.2) 

From this line-element the following set of equations obtains, expressed in terms of the 
extrinsic curvature: 

K t , = —(Ni\ i + N. iu-^p-) 
v b b dt ' 



H = ^g- (KylO* - K 2 + R) = 
Hi = 2^ (KL - K {i ) = 



1 d 
y/gdt 



[y/g {K* - Kg*)] = -N(R l - 



+2N(K ik K J k - K K 3 ) - (N^ 3 - N^g 13 ) + [(K 3 - K g l3 )N l }\i 

-N^K 13 - K g lj ) - N 3 \(K U - K g li ) 



±R g 13 ) - j(K kl K kl - K 2 )tf 3 



(4.4.3a) 
(4.4.3b) 

(4.4.3c) 

(4.4.3d) 
(4.4.3e) 
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This set is equivalent to the ten Einstein's equations. 

In cosmology, we are interested in the class of spatially homogeneous spacetimes, 
characterized by the existence of an m-dimensional isometry group of motions G, acting 
transitively on each surface of simultaneity E t . When m is greater than 3 and there is 
no proper invariant subgroup of dimension 3, the spacetime is of the Kantowski-Sachs 
type [§ and will not concern us further. When m equals the dimension of -which is 
3-, there exist 3 basis one-forms of satisfying: 



^/3 7 °i ° j 



(4.4.4a) 



where Cp are the structure constants of the corresponding isometry group. 

In this case there are local coordinates t, x % such that the line element in ([4.4.2 ) assumes 

the form: 



ds 2 = (iV Q (t)iV a (t) - N 2 {t))dt 2 + 2N a {t)af(x)dx i dt 



(4.4.4b) 
(4.4.4c) 



Latin indices, are spatial with range from 1 to 3. Greek indices, number the different 
basis 1-forms, take values in the same range, and are lowered and raised by j a p, and r ) a/3 
respectively. 

A commitment concerting the topology of the 3-surface, is pertinent here, especially 
in view of the fact that we wish to consider diffeomorphisms [H|; we thus assume that 
G is simply connected and the 3-surface E 4 can be identified with G, by singling out a 
point p of H t , as the identity e, of G. 

If we insert relations ( |4.4.4j ) into equations ( (4.4.3[ ) , we get the following set of ordinary 
differential equations for the Bianchi-Type spatially homogeneous spacetimes: 



E = K% Ki - K 2 + R = 



(4.4.5a) 



K? C. 



<\fi 



(4.4.5b) 



Eg = — NKKq + NRp + 2N P (K« C v 0p - Kl C\ 



where K% = j ap K p p and 



L /3 ^up) 



(4.4.5c) 



P 



2N 



(4.4.6) 



/3k pis 



re A 



(4.4.7) 
(4.4.8) 



When N a = 0, equation ( f4.4.5| c) reduces to the form of the equation given in ||12|| . Equa- 
tion set ( |4.4.5|) , forms what is known as a -complete- perfect ideal; that is, there are no 



15 



integrability conditions obtained from this system. So, with the help of ( [4.4.5 c), (|4.4.6|) , 



( |4.4.7| ), it can explicitly be shown, that the time derivatives of (|4.4.5| a) and ( [4.4. 5[ b) 
vanish identically. The calculation is staightforward -although somewhat lengthy. It 
makes use of the Jacobi identity C^C^ + C^gCp + C^C^ = 0, and its contracted form 

o:.a; o. 

The vanishing of the derivatives of the 4 constrained equations: 
Eq — 0, E a — 0, implies that these equations, are first integrals of equations (|4.4.5| c) 
-moreover, with vanishing integration constants. Indeed, algebraically solving (|4.4.5| a), 
( |4.4.5| b) for N(t),N a (t), respectively and substituting in ( |4.4.5| c) , one finds that in all 
-but Type II and III- Bianchi Types, equations (|4.4.5| c) , can be solved for only 2 of the 6 
accelerations 7 Q/ g present. In Type II and III, the independent accelerations are 3, since 
E a are not independent and thus, can be solved for only 2 of the 3 N a, s. But then in both 
of these linear combination of the N a, s remains arbitrary, and counterbalances 

the extra independent acceleration. Thus, in all Bianchi Types, 4 arbitrary functions of 
time enter the general solution to the set of equations Q4.4.5I). Based on the intuition 



gained from the full theory, one could expect this fact to be a reflection of the only known 
covariance of the theory; i.e. of the freedom to make arbitrary changes of the time and 
space coordinates. 

The rest of this section is devoted to the investigation of the existence, uniqueness, 
and properties of general coordinate transformations -containing 4 arbitrary functions 
of time-, which on the one hand, must preserve the manifest spatial homogeneity, of the 
line element ( f4.4.4| b), and on the other hand, must be symmetries of equations ( |4.4.5| ). 



As far as time reparametrization is concerned the situation is pretty clear: If a transfor- 
mation: 

t -». t = g{t) <=> t = f{t) (4.4.9a) 
is inserted in the line element Q4.4.4J b), it is easily inferred that: 

lapit) - lap{f{t)) EE %{](t) (4.4.9b) 



N(t)^±N(f(t)) d ^ = N(t) 

N a (t) -> N a (f(t))^X^- ee N a {t) 
dt 



(4.4.9c) 
(4.4.9d) 



Accordingly, transforms under (f4.4.9| a) as a scalar and thus ( |4.4.5| a) , ( |4.4.5| b) are 
also scalar equations while ( |4.4.5fc ) gets multiplied by a factor df(t)/dt. Thus, given a 
particular solution to equations ( |4.4.5|) , one can always obtain an equivalent solution, by 
arbitrarily redefining time. Hence, we understand the existence of one arbitrary function 
of time in the general solution to Einstein's equations ([4.4.5 ). In order to understand 
the presence of the rest 3 arbitrary functions of time it is natural to turn our attention 
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to the transformations of the 3 spatial coordinates x l . To begin with, consider the 
transformation: 



i=t& t = t (4.4.10) 
x* = g\x j ,t) x l = f(x j ,t) (4.4.11) 

It is here understood, that our previous assumption concerning the topology of G and 
the identification of E t with G, is valid for all values of the parameter t, for which the 
transformation is to be well defined. 

Under these transformations, the line element ( [4.4.4| b) becomes: 

ds 2 = [(N«N a - N 2 ) + % d -t af{f)a P {f)lapi t) (4.4.12) 

+2of(/)S£ N a (t)]dP (4.4.13) 

+2a? (x)^[N a (t) + aP(x)^ la p(t)}dx m dt (4.4.14) 

+af{x)a P j {x) la p{t)^^ dx m dx n (4.4.15) 

Since our aim, is to retain manifest spatial homogeneity of the line element (|4.4.4| b), we 
have to refer the form of the line element given in ( |4.4. 12| ) to the old basis erf (x) at the 
new spatial point x % . Since erf -both at x % and x 1 -, as well as, dx l jdx\ are invertible 
matrices, there always exists a non-singular matrix A™(x, t) and a triplet P a (x,t), such 
that: 

<^)q^= K&tXW ( 4 - 4 - 16 ) 

af(x)^= P"(x,t) (4.4.17) 

The above relations, must be regarded as definitions, for the matrix Af and the triplet 
P a . With these identifications the line element (f4.4.12| ) assumes the form: 

ds 2 = [(N a N a - N 2 ) + P a (x, t)PP(x, t)7 aP (t) + 2P a (x, t)N a (t)}dt 2 (4.4.18) 
+2Af (x, t)a^(x)[N a (t) + P?(x, i) la(3 (t)]dx m dt (4.4.19) 
+A«(x, t)AP(x, t) laP (i)a^(x)a u n (5;)di m di n (4.4.20) 

If, following the spirit of [TO], we wish the transformation ( ^4.4. 10| ) to be manifest ho- 
mogeneity preserving i.e. to have a well defined, non-trivial action on 7 a /?(t), N(t) and 
N a (t), we must impose the condition that Af(x,t) and P a (x,t) do not depend on the 
spatial point x, i.e. A° = A°(t) and P a = P a (t). Then (|4.4.18|) is written as: 

ds 2 = [(N a N a - N 2 ) + P a P^ aP + 2P a N a ]dt 2 (4.4.21) 
+2A^(x)[N a + P^ a/3 ]dx m dt (4.4.22) 
+A*A^ a/3 a^(x)a»(x)d5; rn dx n => (4.4.23) 

ds 2 = (N a N a - N 2 )dt 2 + 2N a {i)af{x)dx i dt (4.4.24) 
+7a/3(t>f(x)crf(x)rfx i rfF (4.4.25) 
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with the allocations: 



7a/3 = K^l^v 

N a = Ki(N p + P p lpl3 ) and thus N a = S%(N? + P?) 



N = N 



(4.4.26a) 
(4.4.26b) 
(4.4.26c) 



(where S = A -1 ). 

Of course, the demand that A? and P a must not depend on the spatial point x l , changes 
the character of ( 4.4.16Q , from identities, to the following set of differential restrictions 
on the functions defining the transformation: 



df 

dx n ' 



a^f)K a {i)ai{x 



(4.4.27a) 



(4.4.27b) 



Equations ( |4.4.27| ) constitute a set of first-order highly non-linear P.D.E.'s for the un- 
known functions f i . The existence of local solutions to these equations is guaranteed by 
Frobenius theorem |Q as long as the necessary and sufficient conditions: 



d 




)- 9 




\dx m ; 


' dx m 


--( 


'df\ 


d 




sdx m J 


dx m 







hold. Through ( |4.4.27|) and repeated use of (|4.4.4| a), these equations reduce respectively 
to: 



p^r a a v — 1a q 



(4.4.28) 
(4.4.29) 



It is noteworthy that the solutions to ( |4.4.28|) and (|4.4.29| ), -by virtue of ([4.4. 26|) - form 
a group, with composition law: 

= (Ai)?(A a )? 



[p 3 r = (^TJPzf + (p.r 
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where (Ai, Pi) and (A2, P2), are two successive transformations of the form (|4.4.26 ). 



Note also, that a constant automorphism is always a solution of (|4.4.28|) , ( |4.4.29| ); indeed, 
Ajg(t) = A^ and P a (t) = solve these equations. Thus, A^ and P a = can be regarded 
as the remaining gauge symmetry, after one has fully used the arbitrary functions of time, 
appearing in a solution Ajg(i) and P a (t). Consequently one can, at first sight, regard 
all the arbitrary constants encountered when integrating ( |4.4.29| ), as absorbable in the 
shift, since the transformation law for the shift, is then tensorial. This is certainly true, 
as long as there is a non zero initial shift. However, if one has used the independent 
functions of time, in order to set the shift zero, then the constants remaining within 
Al are not absorbable. It is this kind of constants that are explicitly present in "T. 
Christodoulakis, G. Kofinas, E. Korfiatis, G. O. Papadopoulos and A. Paschos, J. Math. 
Phys. 42 (2001) 3580-3608, |gr-qc/0008050| ". Where the solutions to ( gXgg) , (|4.4.29| ) for 



all Bianchi Types are given. A relevant nice discussion, distinguishing between genuine 
gauge symmetries (cf. arbitrary functions of time) and rigid symmetries (cf. arbitrary 
constants), is presented in [TJj]. There a different definition of manifest homogeneity 



preserving diffeomorphisms -stronger than the one adopted in this work- is used, and 
results in only the inner automorphisms being allowed to acquire t dependence. In 
connection to this, it is interesting to observe that (|4.4.28| , |4.4.29|) give essentially the 
same results: notice that 2P tJ- C^ /3 is, by definition, the generator of Inner Automorphisms. 
Thus there is always a \ a p {t) = Exp(2P^C^) G IAut(G) satisfying (|4.4.29|) . If we now 



parameterize the general solution to ( fl.4.28| , (4.4.29[ ) by Ajg(i) = (t)£/|(t) and substitute 



in these relations, we deduce that the matrix U is a constant automorphism. This 
analysis is verified in the explicit solutions to ( jl.4.28|, |4.4.29|) , presented in references 
quoted above. 

4.2 Automorphisms, Invariant Description of 3-spaces, and Quantum 
Cosmology 

As it is well known, the quantum cosmology approximation consists in freezing out all 
but a finite number of degrees of freedom of the gravitational field and quantize the 
rest. This is done by imposing spatial homogeneity. Thus, our -in principle- dynamical 
variables are the scale factors 7^ (t), of some spatially homogeneous geometry. 

The basic object of the theory, is the wave function which must describe the quan- 
tum evolution of the 3-geometry. Consequently, the wave function, will -in principle- 
depend on the 6 7 a /?'s. Hence, a question naturally arises; whether all different ma- 
trices, are characterizing different 3-geometries, or not. The answer to this question, 
involves the A.I.D.s of the previous section, with the difference that time does not con- 
cern us. Thus, the frozen analogue of ( [4.4. 10| ) will lead us to ( |4.4.26| ) (A being now, 



constant) and the integrability condition ( |4.4.28|) . 



So, any two matrices , jfp 5 connected by an element of the automorphism group A^ 

(for an arbitrary albeit given Bianchi Type) i.e. satisfying 7^ = A^A^^L, represent 
the same 3-geometry. 
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The existence of these A.I.D.'s has very important implications for the wave functions 
of a given Bianchi geometry: as we have proven, points in the configuration space - 
spanned by the 7 0( g's, named A- that are related through an automorphism, correspond 
to spatial line elements that are G.C.T. related and thus geometrically identifiable. Thus, 
if we want our wave-functions to depend only on the Geometry of the three-space and 
not on the spatial coordinate system, we must assume them to be annihilated by the 
generators of the entire Automorphism Group and not just by the constraint vector fields 
Hp, which generate only the so-called inner-automorphisms, i.e. we have to demand: 



where: A^ a = (C^V, £^ a ) are the generators of (the connected to the identity component 
of) the Automorphism group and (i) labels the different generators. Depending on the 
particular Bianchi Type, the vector fields (in A) may also include, except of the 
Hp's, the generators of the outer-automorphisms: Ej = £ (j) p la-TQ^— 

Using the method of characteristics, the solutions to the set (|4.4.30| ), can be found 
to have the form: 



where: 



Q l K u , la ,) = —^ (4.4.31a) 



(4.4.31b) 



q\C^ laP ) = — (4.4.31c) 

These three quantities, serve to invariantly describe the geometry and one can prove the 
following relevant preposition (T. Christodoulakis, E. Korfiatis & G.O. Papadopoulos, 
gr-qc/UlU7U50l): 



Let 7„ , 7„ a , G A, and C^ v be the structure constants of a given Bianchi Type. 
If g i (7 (1) ,C) = <f(7( 2 \C) (i = 1,2,3), then there is A such that 7$ = A^A^g? and 
A E Aut(G) i.e. C^A« = A^C« 



It is important to notice that the reduction from a 6-dim configuration space -spanned 
by the 7 aj g- to a space spanned by the g's, is achieved solely by kinematical considerations, 
i.e. the action of G.C.T.'s. Further specification of the theory, may occur only through 
the considerations of the dynamics i.e. the Wheeler-DeWitt equation. 



We close this presentation, by giving two examples. We firstly consider the (see T. 
Christodoulakis, G. O. Papadopoulos, Phys. Lett. B 501 (2001) 264-8): 
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4.2.1 Quantization of the most general Bianchi Type II Vacuum Cosmologies 



In [|IIJ, we had considered the quantization of an action corresponding to the most general 
Bianchi Type II cosmology, i.e. an action giving Einstein's Field Equations, derived from 
the line element: 

ds 2 = (N 2 (t) - N a (t)N a (t))dt 2 + 2N a (t)a?{x)dx i dt + ^ aP (t)a° l (x)a^(x)dx i dx j (4A.32) 
with: 

a a (x) = a° ! (x)dx l 
cr 1 (x) = dx 2 — x 1 dx 3 
a 2 (x) = dx 3 
a 3 (x) = dx 1 
da a (x) = IC^o-P A 
1 



(4.4.33) 



°23 



-r 1 

°32 



sec 
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As is well known [fF7f| , the Hamiltonian is H = N(t)H + N a {t)H a where: 

H = ^L a p^^ + 1 R 
is the quadratic constraint with: 



(4.4.34) 



(4.4.35) 



7 being the determinant of 7 a/ g (the last equality holding only for the Type II case), and: 

H a = (4.4.36) 



are the linear constraints. Note that N appearing in the Hamiltonian, is to be identified 
with N/^fy. 

The quantities H , H a , are weakly vanishing ||, i.e. H « 0, H a w 0. For all class A 
Bianchi Types (C"g = 0), they can be seen to obey the following first-class algebra: 



{H ,H } = 
{H ,H a } = 

{H a , Hp] = - 



(4.4.37) 



which ensures their preservation in time i.e. H « 0, H a « and establishes the 
consistency of the action. 

If we follow Dirac's general proposal || for quantizing this action, we have to turn 
Ho, H a , into operators annihilating the wave function \&. 
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In the Schrodinger representation: 



n «P = _ z a ( 4 - 4 -38) 

satisfying the basic Canonical Commutation Relation (CCR) -corresponding to the clas- 
sical ones: 

= = ^-{SWp + (4.4.39) 

The quantum version of the 2 independent linear constraints has been used to reduce, 
via the method of characteristics |18| , the dimension of the initial configuration space 
from 6 (<y a p) to 4 (combinations of j a/3 ), i.e. * = *(?,7,7i 2 - 7n722, 7i27i3 - 711723), 

where q = C« K C? xla ^r X - 

According to Kuhar's and Hajicek's Jl9[ prescription, the "kinetic" part of Ho is to 



be realized as the conformal Laplacian, corresponding to the reduced metric: 

L a p» ° Z X = g» (4.4.40) 

Ola/3 Ol^u 

where x l , i = 1,2, 3, 4, are the arguments of The solutions had been presented in [|TI| . 



Note that the first-class algebra satisfied by H , H a , ensures that indeed, all components 
of are functions of the af's. The signature of the g l \ is (+,+,—,—) signaling the 
existence of gauge degrees of freedom among the x l, s. 



Indeed, one can prove |20| that the only gauge invariant quantity which, uniquely and 



irreducibly, characterizes a 3-dimensional geometry admitting Type II symmetry group, 



is: 



1 = C a ^ vxla ^l KX (4.4.41) 
An outline of the proof, is as follows: 

Let two hexads 7^ and 7^ be given, such that their corresponding g's, are equal. 



Then according to the result given at the end of the previous section ||20|| , there exists 
an automorphism matrix A (i.e. satisfying C™ A£ = C^KPhl) connecting them, i.e. 



^ l = A^jltJAp. But as it had been shown in the appendix of []2T], this kind of changes 
on j a p, can be seen to be induced by spatial diffeomorphisms. Thus, 3-dimensional Type 
II geometry, is uniquely characterized by some value of q. 



Although for full pure gravity, Kuchaf |22| has shown that there are not other first 



class functions, homogeneous and linear in ir af3 , except H a , imposing the extra symme- 
tries (Type II), allows for such quantities to exist -as it will be shown. We are therefore, 
naturally led to seek the generators of these extra symmetries -which are expected to 
chop off . Such quantities are, generally, called in the literature "Conditional 

Symmetries". 
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The automorphism group for Type II, is described by the following 6 generators -in 
matrix notation and collective form: 



\fo =| k p I (4.4.42) 
with the property: 

c%K = c; a \: + c^x; (4.4.43) 

^From these matrices, we can construct the linear -in momenta- quantities: 

A {1) = X a w l ap n^ (4.4.44) 

Two of these, are the H a ,'s since C?.p correspond to the inner automorphism subgroup 
-designated by the x and y parameters, in A" 7 ^. The rest of them, are the generators of 
the outer automorphisms and are described by the matrices: 




(4.4.45) 



The corresponding -linear in momenta- quantities, are: 

E (I) = e a wla X P (4-4.46) 
The algebra of these -seen as functions on the phase space, spanned by 7^ and tt^-, is: 

{Ei, Ej} = CfjE K 

{E h H a } = -\\P a Hp (4.4.47) 
{E h H } = -2(/t + /i) 7j R 



^From the last of ( |4.4.47|) , we conclude that the subgroup of E^s with the property 



k + /i = 0, i.e. the traceless generators, are first-class quantities; their time derivative 
vanishes. So let: 

E I = {E I : K + fi = 0} (4.4.48) 
Then, the previous statement translates into the form: 

% = => % = a (4.4.49) 

the cj's being arbitrary constants. 

Now, these are -in principle- integrals of motion. Since, as we have earlier seen, 
Ei's along with H a 's, generate automorphisms, it is natural to promote the integrals of 
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motion ( |4.4.49| ), to symmetries -by setting the c/'s zero. The action of the quantum 



version of these Ej's on is taken to be |19 

E f = £ W/3 7a ^ = 1 ^ = V(q, 7) (4.4.50) 

The Wheeler-DeWitt equation now, reads: 

^ TFT - 37 THf + + 5 <?7T - 3 ^7T - 2 ?7* = ( 4 - 4 - 51 ) 

aq^ cry 2 cryag ag c/7 



iVoie that: 



V 2 . = V 2 + -^r-4U = V 2 (4.4.52) 
4(d— 1) 

since we /iawe a 2- dimensional, flat space, with contravariant metric: 

* - ( t ) ^ 

which is Lorentzian. This equation, can be easily solved by separation of variables; 
transforming to new coordinates u = 57 s and v = g 7 , we get the 2 independent equations: 

16u 2 A"(u) + 16uA'(u) - cA(u) =0 , . 

B"(v) + ±B» - (i + £)B(v) = ^ 4 - 4 - 54j 

where c, is the separation constant. Equation fl4.4.51|) , is of hyperbolic type and the 
resulting wave function will still not be square integrable. Besides that, the tracefull 
generators of the outer automorphisms, are left inactive -due to the non vanishing CCR 
with H Q . 

These two facts, lead us to deduce that there must still exist a gauge symmetry, 
corresponding to some -would be, linear in momenta- first-class quantity. Our starting 
point in the pursuit of this, is the third of (|4.4.47| ). It is clear that we need another 
quantity -also linear in momenta- with an analogous property; the trace of tt^, is such 
an object. We thus define the following quantity: 

T = Ej~(K + fi) laf3 7T^ (4.4.55) 

in the phase space -spanned by y a p and n^. It holds that: 

{T,H } = 

{T,H a } = (4.4.56) 
{T,Ej} = 

because of: 



{E / , 7 } = -2( K + /i) 7 
{£/,?} = 

laf3{n a(3 , q} = q 
7^{vr Q/3 ,7} = -37 



(4.4.57) 
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Again -as for E^s-, we see that since T, is first-class, we have that: 

f = T = const = c T (4.4.58) 

another integral of motion. We therefore see, that T has all the necessary properties to 
be used in lieu of the tracefull generator, as a symmetry requirement on In order to 
do that, we ought to set ct zero -exactly as we did with the c/s, corresponding to E^s. 
The quantum version of T, is taken to be: 

f = A ^i" (K+/l)7 "^ (4 ' 4 ' 59) 

Following, Dirac's theory, we require: 

dty dty dty 
™ = A^7„p- {k + n) laP - = (k + yL){q— - = (4.4.60) 

Equation ([4.4.60]) , implies that ty(q,j) = ^(qj) and thus equation (|4.4.51|) , finally, 
reduces to: 

4w 2 V'(w) + 4wV(w) - 2w^ = (4.4.61) 
where, for simplicity, w = qj. The solution to this equation, is: 



* = Ci/ (a/297) + c 2 K (^2q^) (4.4.62) 

where Iq is the modified Bessel function, of the first kind, and K is the modified Bessel 
function, of the second kind, both with zero argument. 

At first sight, it seems that although we have apparently exhausted the symmetries 
of the system, we have not yet been able to obtain a wave function on the space of the 
3-geometries, since $ depends on 57 and not on q only. On the other hand, the fact that 
we have achieved a reduction to one degree of freedom, must somehow imply that the 
wave function found must be a function of the geometry. This puzzle finds its resolution 
as follows. Consider the quantity: 

n = -2 lpa ^ + ±j^L 7 7 "^ CT 4 %p^7«/37 nPa (4 4 63) 

q 

This can also be seen to be first-class, i.e. 

(l = n = const = c n (4.4.64) 

Moreover, it is a linear combination of T, Ei's, and HaS, and thus cq = 0. Now it can 
be verified that fl, is nothing but: 

- + ~) (4-4.65) 



N(ty<y 3q 



25 



So: 

^/g 1 / 3 = $ = constant (4.4.66) 

Without any loss of generality, and since d is not an essential constant of the classical 
system (see [J2B|] and reference [18] therein), we set $ = 1. Therefore: 

¥ = c 1 / (v / 2g 1/3 ) + c 2 K {V2q 1 / 3 ) (4.4.67) 

where Jo is the modified Bessel function, of the first kind, and Kq is the modified Bessel 
function, of the second kind, both with zero argument. 

As for the measure, it is commonly accepted that, there is not a unique solution. A 
natural choice, is to adopt the measure that makes the operator in (|4.4.61| ), hermitian 
-that is: 

li(q) oc q- 1 (4.4.68) 

It is easy to find combinations of C\ and c 2 so that the probability /z(g)|\I/| 2 , be defined. 

Note that putting the constant associated with Q, equal to zero, amounts in restrict- 
ing to a subset of the classical solutions, since cq, is one of the two essential constants 
of Taub's solution. One could keep that constant, at the expense of arriving at a wave 
function with explicit time dependence, since then: 

7 = q~ 1/3 Exp[ J c n N(t)dt] (4.4.69) 

We however, consider more appropriate to set that constant zero, thus arriving at a 
\& depending on q only, and decree its applicability to the entire space of the classical 
solutions. Anyway this is not such a blunder, since \I> is to give weight to all states, 
-being classical ones, or not. 



And the last example (see T. Christodoulakis, G. Gakis & G. O. Papadopoulos, 
gr-qc/0l06065|): 



4.2.2 Conditional Symmetries and the Quantization of Bianchi Type I Vacuum Cos- 
mologies with Cosmological Constant 

Note: The original work, deals with both the cases; the models with a vanishing and those 
with non vanishing cosmological constant. 

The case of Bianchi Type I geometries, has been repeatedly treated in the literature 
-both at the classical level [25] and the quantum level [26]. The main reason, is the 



simplicity brought by the vanishing structure constants, i.e. the high spatial symmetry 
of the model. Thus, the most general of these models, is described by the 6 scale factors 
7a/j(i) and the lapse function N(t) -the shift vector N a (t), being absent due to the non 
existence of the if a 's (linear constraints). The absence of H a s presents -at first sight- the 



26 



complication that no reduction of the initial configuration space, is possible -in contrast 
to what happens in other Bianchi Types p8 |. 



In what follows, we present a complete reduction of the initial configuration space for 
Bianchi Type I geometry, when the cosmological constant is present. A wave function, 
which depends on one degree of freedom, is found. 



As is well known (first of [|17]]) the Hamiltonian of the above system is: 

H = N(t)H + N a (t)H a 

where: 

#o = l -L a p^^ v + 7 A (4.4.70) 

Thus, the only operator which must annihilate the wave function, is H ; and the 
Wheeler-DeWitt equation H Q ^ = 0, will produce a wave function, initially residing on a 
6-dimensional configuration space -spanned by 7 Q/ g's. The discussion however, does not 
end here. If the linear constraints existed, a first reduction of the initial configuration 
space, would take place |T3fl . New variables, instead of the 6 scale factors, would emerge 
-say q l , with i < 6. Then a new "physical" metric would be induced: 

g ij = L^v d J l q * (4.4.71) 

Ola? OJfiv 



According to Kuchaf's and Hajicek's |19| prescription, the "kinetic" part of H would 



have to be realized as the conformal Laplacian (in order for the equation to respect the 
conformal covariance of the classical action), based on the physical metric ( |4.4.71| ). In 
the presence of conditional symmetries, further reduction can take place, a new physical 
metric would then be defined similarly, and the above mentioned prescription, would 
have to be used after the final reduction p8| p2| . 

The case of Bianchi Type I, is an extreme example in which all the linear constraints, 
vanish identically; thus no initial physical metric, exists -another peculiarity reflecting 
the high spatial symmetry of the model under consideration. In compensation, a lot of 
integrals of motion exist and the problem of reduction, finds its solution through the 
notion of "Conditional Symmetries". 

The automorphism algebra of this Type, has been exhaustively treated in the liter- 
ature -see e.g. P7| . The relevant group, is that of the constant, real, 3x3, invertible, 
matrices i.e. GL(3,9ft). The generators of this automorphism group, are (in a collective 
form and matrix notation) the following 9 -one for each parameter: 

/ a (3 5 \ 

>ii)p= [ e ( V ), Je[l,...,9] (4.4.72) 
with the defining property: 

c%K = c; a x: + c^x;. (4.4.73) 
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Exponentiating all these matrices, one obtains the outer automorphism group of Type I. 



For full pure gravity, Kuchaf |22) has shown that there are no other first-class func- 
tions, homogeneous and linear in the momenta, except the linear constraints. If however, 
we impose extra symmetries (e.g. the Bianchi Type I -here considered), such quantities 
may emerge -as it will be shown. We are therefore -according to Dirac ||- justified to 
seek the generators of these extra symmetries, whose quantum-operator analogues will 
be imposed as additional conditions on the wave function. Thus, these symmetries are 
expected to lead us to the final reduction, by revealing the true degrees of freedom. Such 
quantities are, generally, called in the literature "Conditional Symmetries" [22]. From 
matrices ( |4.4.72j ), we can construct the linear -in momenta- quantities: 



E, 



ap 



71 



(4.4.74) 



In order to write analytically these quantities, the following base is chosen: 



Ai 



A i 



A 7 




A, 



A, 



A, 




A, 



A fi 



An 




(4.4.75) 



It is straightforward to calculate the Poisson Brackets of Em with Hq\ 

{E {I)) H Q } = - 1 l\\ a a 

But, it holds that: 

E {I) = {E {Ih H } = - 1 AX a a 



-the last equality emerging by virtue of (4.4.76). Thus 



E, 



{Em, H } = =>- Em = Km = constants, I G [1, 



(4.4.76) 



(4.4.77) 



(4.4.78) 



We therefore conclude that, the first eight quantities Em, are first-class, and thus inte- 
grals of motion. Out of the eight quantities Em, only five are functionally independent 
(i.e. linearly independent, if we allow for the coefficients of the linear combination, to 
be functions of the 7^ s); numerically, they are all independent. 
The algebra of Em can be easily seen to be: 



{E(i),E(j)} 



CjjE(m), 



y,M6[i,...,9] 



(4.4.79) 
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where: 



[A (J) , = CfjX m , I,J,Me [I,... ,9] (4.4.80) 
the square brackets denoting matrix commutation. 

The non vanishing structure constants of the algebra ( [4.4.80 ), are found to be: 



°13 


= 1 


°15 — 




1 


r 7 - 

^16 — 


1 




r 1 - 

°18 - 




1 


°17 — 


-2 


°24 


= 1 


°25 _ 


1 




r<8 

°25 _ 




1 


°26 _ 




1 


u 27 _ 


-1 




= 1 


/-f8 

L/ 34 — 




1 


°35 — 


1 




°37 — 


1 




°38 — 


2 


°46 


= 1 






1 


r~<4 

L/ 48 — 




2 


°57 — 


1 




/"f5 

U 58 — 


-1 


^67 


= 2 


°68 — 


1 





















(4.4.81) 



At this point, in order to achieve the desired reduction, we propose that the quantities 
Em -with / e [1, . . . , 8]- must be promoted to operational conditions acting on the 
requested wave function \l/ -since they are first class quantities and thus integrals of 
motion (see ( [4.4. 78| )). In the Schrodinger representation: 

E {1) * = -zA[ /)a7r/3 ^— = K (J) % I e [1, . . . , 8] (4.4.82) 

In general, systems of equations of this type, must satisfy some consistency conditions 
(i.e. the Frobenius Theorem): 



E {J) ^ = K (J) ^ E {I) E {I) y = K {I) K {J) ^ 



(4.4.83) 



Subtraction of these two and usage of ( f4.4.79| ), results in: 



KfjE {M) y = =>- CfjK^) = (4.4.84) 

i.e. a selection rule for the numerical values of the integrals of motion. Consistency 
conditions ( [4.4.84Q and the Lie Algebra ( |4.4.81| ), impose that Ki = . . . = K 8 = 0. If we 



also had Etg) (as is the case when A = 0) then Kg would remain arbitrary. With this 
outcome, and using the method of characteristics, the system of the five functionally 
independent P.D.E. s ([4.4.82|) , can be integrated. The result is: 

^ = (7) (4.4.85) 

i.e. an arbitrary (but well behaved) function of 7 -the determinant of the scale factor 
matrix. 

A note is pertinent here; from basic abstract algebra, is well known that the basis of 
a linear vector space, is unique -modulo linear mixtures. Thus, although the form of the 



system (U..4.8Q ) is base dependent, its solution ( U-4-8W , is base independent. 

The next step, is to construct the Wheeler-DeWitt equation which is to be solved 
by the wave function (|4.4.85|) . The degree of freedom, is 1; the 5 = 7. According 
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to Kuchaf's proposal ||19|| , upon quantization, the kinetic part of Hamiltonian is to be 
realized as the conformal Beltrami operator - based on the induced physical metric 
-according to (|4.4.71| ), with q = 7: 



g 11 = L^p-p- = L^VY* = -3 7 2 (4-4.86) 



In the Schrodinger representation: 



\L a ,^^ -> -\u\ (4.4.87) 



where: 



ni = n'=—= d^y/gTi 9 LL d,} (4.4.88) 

is the 1-dimensional Laplacian based on gu {g ll gu = 1). Note that in 1-dimension the 
conformal group is totally contained in the G.C.T. group, in the sense that any conformal 
transformation of the metric can not produce any change in the -trivial- geometry and 
is thus reachable by some G.C.T. Therefore, no extra term in needed in (J4.4.88), as it 



can also formally be seen by taking the limit d — 1, R = in the general definition: 

(d-2) 

A(d-iy 



D 2 = D 2 + id 2) R= g2 



Thus: 



H ^H = -i(-3 7 2 |^ - 37^) + A7 (4.4.89) 



So, the Wheeler-DeWitt equation -by virtue of ( 14.4.851) -, reads: 

2 



H ty = 7 V + 7*' + - 7 A^ = (4.4.90) 

3 



The general solution to this equation, is: 



*( 7 ) = Cl J (2y ^) + c 2 F (2^/^) (4.4.91) 

where J n and Y n , are the Bessel Functions of the first and second kind respectively -both 
with zero argument- and c±,C2, arbitrary constants. 

An important element for selecting the measure, is the conformal covariance; the 
supermetric L a/3fJ-l/ is known only up to rescalings, because instead of N(t) one can take 
any N(t) = N(t)e~ 2uj (with cu = uj(j aP )) and consequently T^ 1 \t) = L a ^ u (t)e 2u . 
This property, is also inherited to the physical metric ( |4.4.71| ) and is the reason for the 
Kuchaf's recipe, adopted in this work. 
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It is natural that the proposed measure /i, must be such that the probability den- 
sity (j, | $ | 2 , be invariant under these scalings. Recalling that ^ = ^/e*- 2- x> ^/ 2 , 
we conclude that fi must scale as /Z = fie^ -2 ^ . The natural measure under which 
the Wheeler-DeWitt operator is hermitian, is ^/Det (physical metric), but it scales as 



Det (physical metric) = y/Det (physical metric)e x>w . 
We are thus after a quantity £ -preferably constant, (so that the hermiticity is 
preserved)- which scales as £ = £e~ 2w . It is not difficult to imagine such a quantity: 
The inverse of any product of E^ a p with E(j^ u (where E^ a p = l/2(A*7 Kj g + (a <-> j3))) 
has the desired property. Indeed the E^ s do not scale at all, while the supermetric 
scales as mentioned before. The group metric Qjj = Cf s Cj F can serve to close the group 
indices of E^ a p. So, we arrive at the quantity: 



L^e IJ E (n „ B E 



(4.4.92) 



(where G /J is the inverse of the group metric) having the desired property and being 
also a constant. Using the Lie algebra ( |4.4.81| ), one obtains: 



Br.7 = CfaCi 



J IS^JF 
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(4.4.93) 



Thus: 
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(4.4.94) 



After a straightforward calculation, one finds that: 



(4.4.95) 



The product of £ with the respective natural measure, defines the final expression for 
the measure /i. 
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It is fair to say that the problem of selection of the "correct" measure, is not yet 
resolved; it is a reflection of the problem of time in Quantum Gravity |29| . 

Another issue that has not been touched upon is the problem of selecting a unique 
wave faction. In the path integral approach, to quantum cosmology, there is the Hartle- 
Hawking "no boundary proposal" f30fl . In the canonical approach, there are various 
forms of the Vilenkin proposal |31] . 



Finally, the problem of decoherence (i.e. of reconstruction of classical trajectories, 
from the knowledge of the wave function), has occupied several workers in the field [B3 . 
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